Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 15 August 2011 (MN WT^ style file v2.2) 



The WiggleZ Dark Energy Survey: mapping the 
distance-redshift relation with baryon acoustic oscillations 

Chris Blake^*, Eyal A. Kazin^, Florian Beutler^, Tamara M. Davis^'^, 
David Parkinson^, Sarah Brough^, Matthew Colless^, Carlos Contreras^, 
Warrick Couch^, Scott Croom^, Darren Croton\ Michael J. Drinkwater^, 
Karl Forster^, David Gilbank^, Mike Gladders^*^, Karl Glazebrook^, 
Ben Jelliffe^, Russell J. Jurek^^, I-hui Li^, Barry Madore^^, 

D. Christopher Martin^, Kevin Pimbblet^^, Gregory B. Poole^, Michael Pracy^'^^, 
Rob Sharp6'i4, Emily Wisnioski\ David Woods^^ Ted K. Wyder^ and H.K.C. Yee^^ 

^ Centre for Astrophysics & Supercomputing, Swinburne University of Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia 

^ Center for Cosmology and Particle Physics, New York University, 4 Washington Place, New York, NY 10003, United States 

^ International Centre for Radio Astronomy Research, University of Western Australia, 35 Stirling Highway, Perth WA 6009, Australia 

^ School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 

^ Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark 
^ Australian Astronomical Observatory, P.O. Box 296, Epping, NSW 1710, Australia 

Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 
* California Institute of Technology, MC 278-17, 1200 East California Boulevard, Pasadena, CA 91125, United States 
^ Astrophysics and Gravitation Group, Department of Physics and Astronomy, University of Waterloo, Waterloo, ON N2L 3C1, Canada 

Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, United States 

Australia Telescope National Facility, CSIRO, Epping, NSW 1710, Australia 
^'^ Observatories of the Carnegie Institute of Washington, 813 Santa Barbara St., Pasadena, CA 91101, United States 

School of Physics, Monash University, Clayton, VIC 3800, Australia 

Research School of Astronomy & Astrophysics, Australian National University, Weston Creek, ACT 2611, Australia 
Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T IZl, Canada 
^® Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada 



15 August 2011 



ABSTRACT 

We present measurements of the baryon acoustic peak at redshifts z — 0.44, 0.6 and 
0.73 in the galaxy correlation function of the final dataset of the WiggleZ Dark Energy 
Survey. We combine our correlation function with lower-redshift measurements from 
the 6-degree Field Galaxy Survey and Sloan Digital Sky Survey, producing a stacked 
survey correlation function in which the statistical significance of the detection of 
the baryon acoustic peak is 4.9-0' relative to a zero-baryon model with no peak. We 
fit cosmological models to this combined baryon acoustic oscillation (BAO) dataset 
comprising six distance-redshift data points, and compare the results to similar fits to 
the latest compilation of supernovae (SNe) and Cosmic Microwave Background (CMB) 
data. The BAO and SNe datasets produce consistent measurements of the equation-of- 
state w of dark energy, when separately combined with the CMB, providing a powerful 
check for systematic errors in either of these distance probes. Combining all datasets 
we determine w ~ —1.03 ± 0.08 for a flat Universe, consistent with a cosmological 
constant model. Assuming dark energy is a cosmological constant and varying the 
spatial curvature, we find ilk — — 0.004 ± 0.006. 

Key words: surveys, large-scale structure of Universe, cosmological parameters, 
distance scale, dark energy 
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Measurements of the cosmic distance-redshift relation have 
always constituted one of the most important probes of the 
cosmological model. Eighty years ago such observations pro- 
vided evidence that the Universe is expanding; more recently 
they have convincingly suggested that this expansion rate 
is accelerating. The distance-redshift relation depends on 
the expansion history of the Universe, which is in turn gov- 
erned by its physical contents including the properties of 
the "dark energy" which has been hypothesized to be driv- 
ing the accelerating expansion. One of the most important 
challenges in contemporary cosmology is to distinguish be- 
tween the different possible physical models for dark energy, 
which include a material or scalar field smoothly filling the 
Universe with a negative equation-of-statc, a modification 
to the laws of gravity at large cosnuc scales, or the effects of 
inhomogeneity on cosmological observations. Cosmological 
distance measurements provide one of the crucial observa- 
tional datasets to help distinguish between these different 
models. 

One of the most powerful tools for mapping the 
distance-redshift relation is Type la supernovae (SNe la). 
About a decade ago, observations of nearby and distant 
SNe la provided some of the most compelling evidence that 
the expansion rate of the Universe is accelerating (Riess et 
al. 1998, Perlmutter et al. 1999), in agreement with ear- 
lier suggestions based on comparisons of the Cosmic Mi- 
crowave Background (CMB) and large-scale structure data 
(Efstathiou, Sutherland & Maddox 1990, Krauss & Turner 
1995, Ostriker & Steinhardt 1995). Since then the sample of 
SNe la available for cosmological analysis has grown impres- 
sively due to a series of large observational projects which 
has populated the Hubble diagram across a range of red- 
shifts. These projects include the Nearby Supernova Fac- 
tory (Copin et al. 2006), the Center for Astrophysics SN 
group (Hicken et al. 2009) , the Carnegie Supernova Project 
(Hamuy et al. 2006) and the Palomar Transient Factory 
(Law et al. 2009) at low redshifts z < 0.1; the Sloan Dig- 
ital Sky Survey (SDSS) supernova survey (Kessler et al. 
2009) at low-to- intermediate redshifts 0.1 < z < 0.3; the Su- 
pernova Legacy Survey (Astier et al. 2006) and ESSENCE 
(Wood-Vasey et al. 2007) projects at intermediate redshifts 
0.3 < z < 1.0; and observations by the Hubble Space Tele- 
scope at high redshifts z > 1 (Riess et al. 2004, 2007; Daw- 
son et al. 2009). These supernovae data have been collected 
and analyzed in a homogeneous fashion in the "Union" SNe 
compilations, initially by Kowalski et al. (2008) and most re- 
cently by AmanuUah et al. (2010) in the "Union 2" sample 
of 557 SNe la. 

The utility of these supernovae datasets is now lim- 
ited by known (and potentially unknown) systematic errors 
which could bias cosmological fits if not handled correctly. 
These systematics include redshift-dependent astrophysical 
effects, such as potential drifts with redshift in the relations 
between colour, luminosity and light curve shape owing to 
evolving SNe la populations, and systematics in analysis 
such as the fitting of light curves, photometric zero-points, 
K-corrections and Malmquist bias. Although these system- 
atics have been treated very thoroughly in recent supernovae 
analyses, it is clearly desirable to cross-check the cosmolog- 
ical conclusions with other probes of the distance-redshift 
relation. 

A very promising and complementary method for map- 



ping the distance-redshift relation is the measurement of 
baryon acoustic oscillations (BAOs) in the large-scale clus- 
tering pattern of galaxies, and their application as a cos- 
mological standard ruler (Eisenstein, Hu & Tegmark 1998, 
Cooray et al. 2001, Eisenstein 2003, Blake & Glazebrook 
2003, Seo & Eisenstein 2003, Linder 2003, Hu & Haiman 
2003). BAOs correspond to a preferred length scale im- 
printed in the distribution of photons and baryons by the 
propagation of sound waves in the relativistic plasma of 
the early Universe (Peebles & Yu 1970, Sunyaev & Zel- 
dovitch 1970, Bond & Efstathiou 1984, Holtzman 1989, Hu 
& Sugiyama 1996, Eisenstein & Hu 1998). This length scale, 
which corresponds to the sound horizon at the baryon drag 
epoch denoted by rs{zd), may be predicted very accurately 
by measurements of the CMB which yield the physical mat- 
ter and baryon densities that control the sound speed, ex- 
pansion rate and recombination time in the early Universe: 
the latest determination is rs{zd) = 153.3 ± 2.0 Mpc (Ko- 
matsu et al. 2009). In the pattern of late-time galaxy clus- 
tering, BAOs manifest themselves as a small preference for 
pairs of galaxies to be separated by rs{zd), causing a distinc- 
tive "baryon acoustic peak" to be imprinted in the 2-point 
galaxy correlation function. The corresponding signature in 
Fourier space is a series of decaying oscillations or "wiggles" 
in the galaxy power spectrum. 

Measurement of BAOs has become an important moti- 
vation for galaxy redshift surveys in recent years. The small 
amplitude of the baryon acoustic peak, and the large size 
of the relevant scales, implies that cosmic volumes of order 
1 Gpc^ must be mapped with of order 10* galaxies to en- 
sure a robust detection (Tegmark 1997, Blake & Glazebrook 
2003, Glazebrook & Blake 2005, Blake et al. 2006). Signif- 
icant detections of BAOs have now been reported by three 
independent galaxy surveys, spanning a range of redshifts 
z < 0.6: the SDSS, the WiggleZ Dark Energy Survey, and 
the 6-degree Field Galaxy Survey (6dFGS). 

The most accurate BAO measurements have been ob- 
tained by analyzing the SDSS, particularly the Luminous 
Red Galaxy (LRG) component. Eisenstein et al. (2005) re- 
ported a convincing detection of the acoustic peak in the 
2-point correlation function of the SDSS Third Data Re- 
lease (DR3) LRG sample with effective redshift z = 0.35. 
Percival et al. (2010) performed a power-spectrum analy- 
sis of the SDSS DR7 dataset, considering both the main 
and LRG samples, and measured the distance-redshift re- 
lation at both z = 0.2 and z = 0.35 with ~ 3% accuracy 
in units of the standard ruler scale. Other studies of the 
SDSS LRG sample, producing broadly similar conclusions, 
have been undertaken by Hutsi (2006) , Percival et al. (2007) , 
Sanchez et al. (2009) and Kazin et al. (2010a). These studies 
of SDSS galaxy samples built on hints of BAOs reported by 
the 2-degree Field Galaxy Redshift Survey (Percival et al. 
2001, Cole et al. 2005) and combinations of smaller datasets 
(Miller et al. 2001). There have also been potential BAO de- 
tections in photometric-redshift catalogues from the SDSS 
(Blake et al. 2007, Padmanabhan et al. 2007, Crocce et al. 
2011), although the statistical significance of these measure- 
ments currently remains much lower than that which can be 
obtained using spectroscopic redshift catalogues. 

These BAO detections have recently been supplemented 
by new measurements from two different surveys, which 
have extended the redshift coverage of the standard-ruler 
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technique. In the low-rodshift Universe the 6dFGS has re- 
ported a BAO detection at 2 = 0.1 (Beutler et al. 2011). 
This study produced a ~ 5% measurement of the standard 
ruler scale and a new determination of the Hubble constant 
Hq. At higher rcdshifts the WigglcZ Survey has quantified 
BAOs at z = 0.6, producing a ~ 4% measurement of the 
baryon acoustic scale (Blake et al. 2011). Taken together, 
these different galaxy surveys have demonstrated that BAO 
standard-ruler measurements are self-consistent with the 
standard cosmological model established from CMB obser- 
vations, and have yielded new, tighter constraints on cosmo- 
logical parameters. 

The accuracy with which BAOs may be used to deter- 
mine the distancc-rcdshift relation using current surveys is 
limited by statistical rather than systematic errors (in con- 
trast to observations of SNe la). The measurement error in 
the large-scale correlation function, which governs how accu- 
rately the preferred scale may be extracted, is determined by 
the volume of the large-scale structure mapped and the num- 
ber density and bias of the galaxy tracers. There are indeed 
potential systematic errors associated with fitting models to 
the BAO signature, which are caused by the modulation of 
the pattern of linear clustering laid down in the high-redshift 
Universe by the non-linear scale-dependent growth of struc- 
ture, the distortions apparent when the signal is observed 
in redshift-space and the bias with which galaxies trace the 
network of matter fluctuations. However, the fact that the 
BAOs arc imprinted on large, linear and quasi-linear scales 
of the clustering pattern means that these non-linear, sys- 
tematic distortions are amenable to analytical or numerical 
modelling and the leading-order effects arc well-understood 
(Eisenstein, Seo & White 2007, Croccc & Scoccimarro 2008, 
Matsubara 2008, Sanchez, Baugh & Angulo 2008, Smith, 
Scoccimarro & Sheth 2008, Seo et al. 2008, Padmanabhan 
& White 2009). As such, BAOs in current datasets are be- 
lieved to provide a robust probe of the cosmological model, 
relatively free of systematic error and dominated by statisti- 
cal errors. In this sense they provide a powerful cross-check 
of the distance-redshift relation mapped by supernovae. 

In this study we report our final analysis of the baryon 
acoustic peak from the angle-averaged correlation function 
of the completed WiggleZ Survey dataset, in which we 
present distance-scale measurements as a function of redshift 
between z = 0.44 and z = 0.73, including a covariance ma^ 
trix which may be applied in cosmological parameter fits. We 
also present a new measurement of the correlation function 
of the SDSS-LRG sample. We stack the 6dFGS, SDSS-LRG 
and WiggleZ correlation functions to produce the highest- 
significance detection to date of the baryon acoustic peak 
in the galaxy clustering pattern. We perform cosmological 
parameter fits to this latest BAO distance dataset, now com- 
prising data points at six different redshifts. By comparing 
these fits with those performed on the latest compilation 
of SNe la, we search for systematic disagreements between 
these two important probes of the distance-redshift relation. 

The structure of our paper is as follows: in Section 2 
we summarize the three galaxy spectroscopic redshift sur- 
vey datasets which have provided the most significant BAO 
measurements. In Section 3 we outline the modelling of the 
baryon acoustic peak applied in this study. In Section 4 we 
report the measurement and analysis of the final WiggleZ 
Survey correlation functions in redshift slices, and in Sec- 



tion 5 we present the new determination of the correlation 
function of SDSS LRGs. In Section 6 we construct a stacked 
galaxy correlation function from these surveys and analyze 
the statistical significance of the BAO detection contained 
therein. In Section 7 we perform cosmological parameter fits 
to various combinations of BAO, SNe la and CMB data, and 
we list our conclusions in Section 8. 



2 DATASETS 

2.1 The WiggleZ Dark Energy Survey 

The WiggleZ Dark Energy Survey (Drinkwater et al. 2010) 
is a large-scale galaxy redshift survey of bright emission- 
line galaxies which was carried out at the Anglo- Australian 
Telescope between August 2006 and January 2011 using the 
AAOmega spectrograph (Saunders et al. 2004, Sharp et al. 
2006). Targets were selected via joint ultraviolet and opti- 
cal magnitude and colour cuts using input imaging from the 
Galaxy Evolution Explorer (GALEX) satellite (Martin et al. 
2005), the Sloan Digital Sky Survey (SDSS; York et al. 2000) 
and the 2nd Red Cluster Sequence (RCS2) Survey (Gilbank 
et al. 2011). The survey is now complete, comprising of order 
200,000 rcdshifts and covering of order 800 deg^ of equato- 
rial sky. In this study we analyzed a galaxy sample drawn 
from our final set of observations, after cuts to maximize the 
contiguity of each survey region. The sample includes a total 
oi N = 158,741 galaxies in the redshift range 0.2 < 2; < 1.0. 



2.2 The 6-degree Field Galaxy Survey 

The 6-dcgrcc Field Galaxy Survey (6dFGS, Jones et al. 
2009) is a combined redshift and peculiar velocity survey 
covering nearly the entire southern sky with the exception 
of a 10° band along the Galactic plane. Observed galaxies 
were selected from the 2MASS Extended Source Catalog 
(Jarrett et al. 2000) and the redshifts were obtained with 
the 6-degrcc Field (6dF) multi-fibre instrument at the U.K. 
Schmidt Telescope between 2001 and 2006. The final 6dFGS 
sample contains 75,117 galaxies distributed over ~ 17,000 
deg^ with a mean redshift of 2 = 0.052. The analysis of 
the baryon acoustic peak in the 6dFGS (Beutler et al. 2011) 
utilized all galaxies selected to K < 12.9. We provide a sum- 
mary of this BAO measurement in Section 6.1. 



2.3 The Sloan Digital Sky Survey Luminous Red 
Galaxy sample 

The SDSS included the largest-volume spectroscopic LRG 
survey to date (Eisenstein et al. 2001). The LRGs were 
selected from the photometric component of SDSS, which 
imaged the sky at high Galactic latitude in five passbands 
u, g, r, i and 2 (Fukugita et al. 1996, Gunn et al. 1998) us- 
ing a 2.5m telescope (Gunn et al. 2006). The images were 
processed (Lupton et al. 2001, Stoughton et al. 2002, Pier 
et al. 2003, Ivezic et al. 2004) and calibrated (Hogg et al. 
2001, Smith et al. 2002, Tucker et al. 2006), allowing se- 
lection of galaxies, quasars (Richards et al. 2002) and stars 
for follow-up spectroscopy (Eisenstein et al. 2001, Strauss et 
al. 2002) with twin fibre-fed double spectographs. Targets 
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were assigned to plug plates according to a tiling algorithm 
ensuring nearly complete samples (Blanton et al. 2003). 

The LRG sample serves as a good tracer of matter be- 
cause these galaxies are associated with massive dark mat- 
ter halos. The high luminosity of LRGs enables a large vol- 
ume to be efficiently mapped, and their spectral uniformity 
makes them relatively easy to identify. In this study we ana^ 
lyze similar LRG catalogues to those presented by Kazin et 
al. (2010a, 2010b)t, to which we refer the reader for full de- 
tails of selection and systcmatics. In particular, in this study 
we focus on the sample DR7-Full, which corresponds to all 
LRGs in the redshift range 0.16 < a; < 0.44 and absolute 
magnitude range —23.2 < Mg < —21.2. The sky coverage 
and redshift distributions of the LRG samples are presented 
in Figures 1 and 2 of Kazin et al. (2010a). DR7-Full in- 
cludes 89,791 LRGs with average redshift (z) = 0.314, cover- 
ing total volume 1.2 Gpc^ with average number density 
8 X 10-^ Mpc"^. 



3 MODELLING THE BARYON ACOUSTIC 
PEAK 

In this Section we summarize the two models we fitted to 
the new baryon acoustic peak measurements presented in 
this study. These models describe the quasi-linear effects 
which cause the acoustic feature and correlation function 
shape to deviate from the linear-theory prediction. There 
are two main aspects to model: a damping of the acous- 
tic peak caused by the displacement of matter duo to bulk 
flows, and a distortion in the overall shape of the cluster- 
ing pattern due to the scale-dependent growth of structure 
(Eisenstein et al. 2007, Croccc & Scoccimarro 2008, Matsub- 
ara 2008, Sanchez et al. 2008, Smith et al. 2008, Seo et al. 
2008, Padmanabhan & White 2009). Our models are char- 
acterized by four variable parameters: the physical matter 
density Qmh'^ (where Qm is the matter density relative to the 
critical density and h — _ffo/[100 kms~^ Mpc~^] is the Hub- 
ble parameter), a scale distortion parameter a, a physical 
damping scale cr„, and a normalization factor b^. The mod- 
els for the correlation function ^modei in terms of separation 
s can be written in the form 

^model(s) = ^Rd{^inh'^,(Tv,as). (1) 

The physical matter density Clmh^ determines (to first order) 
both the overall shape of the matter correlation function and 
the length scale of the standard ruler, by determining the 
physics before recombination. The scale distortion param- 
eter a relates the distance-redshift relation at the effective 
redshift of the sample to the fiducial value used to construct 
the correlation function measurement, in terms of the Dv 
parameter (Eisenstein et al. 2005, Padmanabhan & White 
2008, Kazin, Sanchez & Blanton 2011): 



Dv{Zes) = aDv,fid{Zes), 



(2) 



where Dv is a composite of the physical angular-diameter 
distance Da{z) and Hubble parameter H{z), which respec- 
tively govern tangential and radial separations in a cosmo- 
logical model: 

t These catalogues and the associated survey mask are publicly 
available at http://cosmo.nyu.edu/~eak306/SDSS-LRG.html 



Dv{z) = 



{l+zfDA{zf 



H{z) 



1/3 



(3) 



The damping scale <t„ quantifies the typical displacement of 
galaxies from their initial locations in the density field due to 
bulk flows, resulting in a "washing-out" of the baryon oscilla- 
tions at low redshift. The normalization factor , marginal- 
ized in our analysis, models the effects of linear galaxy bias 
and large-scale redshift-space distortions. 



3.1 Default correlation function model 

In our first, default, model we constructed the fiducial cor- 
relation function ^fld in Equation 1 in a similar manner to 
Eisenstein et al. (2005) and Blake et al. (2011). First, we 
generated a linear power spectrum Pi^{k) as a function of 
wavenumber k for a given fim/i^ using the CAMB software 
package (Lewis, Challinor & Lasenby 2000). We fixed the 
values of the other cosmological parameters using a fidu- 
cial model consistent with the latest fits to the Cosmic Mi- 
crowave Background (Komatsu et al. 2011): Hubble param- 
eter h = 0.71, physical baryon density Q,hh^ = 0.0226, 
primordial spectral index rig = 0.96 and normalization 
£78 = 0.8. We also used the fitting formulae of Eisenstein 
& Hu (1998) to generate a corresponding "no- wiggles" ref- 
erence spectrum Pref (fe), possessing a similar shape to Pi^{k) 
but with the baryon oscillation component deleted, which we 
also use in the clustering model as explained below. 

We then incorporated the damping of the baryon acous- 
tic peak caused by the displacement of matter due to bulk 
flows (Eisenstein et al. 2007, Crocce & Scoccimarro 2008, 
Matsubara 2008) by interpolating between the linear and 
reference power spectra using a Gaussian damping term 
g{k) = cxp {-k^aly. 



Pdamped(fc) = g{k) Pl (fc) + [1 - g{k)] P,,i{k) . 



(4) 



The magnitude of the damping coefHcient cr„ can be esti- 
mated for a given value of fim/i^ using the first-order pre- 
diction of perturbation theory (Crocce & Scoccimarro 2008): 



2 



1 



I P^{k)dk. 



(5) 



However, this relation provides only an approximation to the 
true non- linear damping (Taruya et al. 2010), and we chose 
to marginalize over cr„ as a free parameter in our analysis. 
We note that it„ is closely related to the parameter fc* de- 
fined by Sanchez et al. (2008), in the sense that al — l/2k^. 

We included the boost in small-scale clustering power 
due to the non-linear scale-dependent growth of structure 
using the "halofit" prescription of Smith et al. (2003), as 
applied to the no-wiggles reference spectrum: 



i^L(fc) ^ 



Pref , halofit 

(fc) 

Pref (A;) 



^ ■fciamped(^)- 



(6) 



Finally, we transformed P^i^ik) into the correlation function 
appearing in Equation 1: 



^fid(s) 



dkk'^Ptn^ik) 



sin (ks) 



ks 



(7) 
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3.2 Comparison correlation function model 

The second, comparison model we considered for the fidu- 
cial correlation function ^fid was motivated by perturbation 
theory (Crocce & Scoccimarro 2008, Sanchez at al. 2008): 



^fid(s) = ^l(s) (8) exp (-s^/2o-^) + Amc 



ds 



ii{s)- (8) 



In this relation 5l(s) is the linear model correlation func- 
tion corresponding to the linear power spectrum Ph{k). The 
symbol (g) denotes convolution by the Gaussian damping cr^ , 
which we evaluated as 



^l(s) ® exp {~s^/2al) 



dkk^ Phik) exp(— fc^cr^) 



sin (ks) 



ks 



(9) 



and ^1 is defined by Equation 32 in Crocce & Scoccimarro 
(2008): 



6(«) 



dkkPi,{k)ji{ks), 



(10) 



where ji{x) is the spherical Bessel function of order 1. 
^Mc = 1 (fixed in our analysis) is a "mode-coupling" term 
that restores the small-scale shape of the correlation func- 
tion and causes a slight shift in the peak position compared 
to the linear-theory prediction. The model of Equation 8 
has been shown to yield unbiased results in baryon acoustic 
peak fits by Sanchez et al. (2008, 2009). 



4 WIGGLEZ BARYON ACOUSTIC PEAK 
MEASUREMENTS IN REDSHIFT SLICES 

In this Section we describe our measurement and fitting 
of the baryon acoustic peak in the WiggleZ Survey galaxy 
correlation function in three overlapping redshift ranges: 
0.2 < z < 0.6, OA < z < 0.8 and 0.6 < 2 < 1.0. Our 
methodology closely follows that employed by Blake et al. 
(2011), to which we refer the reader for full details. 

4.1 Correlation function mcEisurements 

We measured the angle-averaged 2-point correlation func- 
tion ^(s) for each WiggleZ survey region using the Landy- 
Szalay (1993) estimator: 



DD{s) - 2 DR(s) + RR(s) 
RR{s) 



(11) 



where DD{s), DR(s) and RR{s) arc the data-data, data- 
random and random-random weighted pair counts in sep- 
aration bin s, where each random catalogue contains the 
same number of galaxies as the real dataset. We assumed a 
fiducial fiat ACDM cosmological model with matter density 
ilm ~ 0.27 to convert the galaxy redshifts and angular posi- 
tions to spatial co-moving co-ordinates. In the construction 
of the pair counts each data or random galaxy i was assigned 
a weight Wi = 1/(1 + njPo), where rii is the survey number 
density at the location of the ith galaxy (determined by av- 
eraging over many random catalogues) and Po = 5000 h^"^ 
Mpc^ is a characteristic power spectrum amplitude at the 
physical scales of interest. The DR and RR pair counts were 
determined by averaging over 10 random catalogues, which 



were constructed using the selection-function methodology 
described by Blake et al. (2010). We measured the correlar 
tion function in 10 Mpc separation bins in three over- 
lapping redshift slices 0.2 < z < 0.6, 0.4 < z < 0.8 and 
0.6 < 2 < 1.0. The effective redshift Zcs of the correlation 
function measurement in each slice was determined as the 
weighted mean redshift of the galaxy pairs in the separation 
bin 100 < s < 110 Mpc, where the redshift of a pair is 
simply the average (21 + 22)/2, and the weighting is ^1^2 
where Wi is defined above. For the three redshift slices in 
question we obtained values ZeB = 0.44, 0.60 and 0.73. 

We determined the covariance matrix of the correla- 
tion function measurement in each survey region using an 
ensemble of 400 lognornial realizations, using the method 
described by Blake et al. (2011). Lognormal realizations 
provide a reasonably accurate galaxy clustering model for 
the linear and quasi-linear scales which are important for 
the modelling of baryon oscillations. They arc more reliable 
than jack-knife errors, which provide a poor approximation 
for the correlation function variance on BAO scales because 
the pair separations of interest are usually comparable to 
the size of the jack-knife regions, which are then not strictly 
independent. We note that the lognormal covariance ma- 
trix only includes the effects of the survey window function, 
and neglects the covariance due to the non-linear growth 
of structure and redshift-space effects. The full non-linear 
covariance matrix may be studied with the aid of a large 
set of N-body simulations (Rimes & Hamilton 2005, Taka- 
hashi et al. 2011). Work is in progress to construct such 
a simulation set for WiggleZ galaxies, although this is a 
challenging computational problem because the typical dark 
matter haloes hosting the star-forming galaxies mapped by 
WiggleZ arc ~ 20 times lower in mass than the Luminous 
Red Galaxy sample described in Section 5, requiring high- 
resolution large-volume simulations. However, we note that 
Takahashi et al. (2011) demonstrated that the impact of us- 
ing the full non-linear covariance matrix on the accuracy 
of extraction of baryon acoustic oscillations is small, so we 
do not expect our measurements to be compromised signifi- 
cantly through using lognormal realizations to estimate the 
covariance matrix. 

We combined the correlation function measurements 
and corresponding covariance matrices for the different sur- 
vey regions using optimal inverse-variance weighting in each 
separation bin (see equations 8 and 9 in White et al. 2011): 



■comb = 



comb 



E or' 



(12) 



(13) 



In these equations, and ^^^^^ are vectors represent- 
ing the correlation function measurements in region n and 
the optimally-combined correlation function, and C_ and 

n 

C_ are the covariance matrices corresponding to these 

comb 

two measurement vectors (with inverses C7^^ and ). 

n comb 

This method produces an almost identical result to combin- 
ing the individual pair counts and then estimating the corre- 
lation function using Equation 11. The combined correlation 
functions in the three redshift slices are displayed in Figure 
1, together with a total WiggleZ correlation function for the 
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Figure 1. Measurements of the galaxy correlation function ^(s), combining different WiggleZ survey regions, for the redshift ranges 
0.2 < 2 < 0.6, 0.4 < z < 0.8, 0.6 < z < 1.0 and 0.2 < z < 1.0, plotted in the combination ^(s) where s is the co-moving redshift-space 
separation. The best-fitting clustering models in each case, varying the parameters flmh'^, ct, and h"^ as described in Section 3, are 
overplotted as the solid lines. Significant detections of the baryon acoustic peak are obtained in each separate redshift slice. 



whole redshift range 0.2 < z < 1.0 which was constructed 
by combining the separate measurements for 0.2 < z < 0.6 
and 0.6 < z < 1.0. The corresponding lognormal covariance 
matrices for each measurement are shown in Figure 2. 

4.2 Parameter fits 

We fitted the first, defauh correlation function model de- 
scribed in Section 3 to the WiggleZ measurements in red- 
shift slices 0.2 < 2 < 0.6, 0.4 < z < 0.8 and 0.6 < z < 1.0, 
varying fim/i'^, Q, cr„ and }? . Our default fitting range was 
10 < s < 180 /i"^ Mpc (following Eisenstein et al. 2005), 
where 10 Mpc is an estimate of the minimum scale of 
validity for the quasi- linear theory described in Section 3. 
This minimum scale is a quantity which depends on the sur- 
vey redshift and galaxy bias (which control the amplitude of 
the non-linear, scale-dependent contributions to the shape of 
the correlation function) together with the signal-to-noise of 
the measurement. When fitting Equation 7 to the WiggleZ 
Survey correlation function we find no evidence for a system- 
atic variation in the derived BAO parameters when we vary 
the minimum fitted scale over the range 10 < Smin < 50 h"^ 
Mpc. 

We minimized the statistic using the full data co- 
variance matrix derived from lognormal realizations. The 
fitting results, including the marginalized parameter mea- 
surements, are displayed in Table 1. The minimum values 
of for the model fits in the three redshift slices were 
11.4, 10.1 and 13.7 for 13 degrees of freedom, indicating 



that our model provides a good fit to the data. The best- 
fitting scale distortion parameters, which provide the value 
of Dvizafi) for each redshift slice, are all consistent with the 
fiducial distance-redshift model (a flat ACDM Universe with 
fJm = 0.27) with marginalized errors of 9.1%, 6.5% and 6.4% 
in the three redshift slices. The best-fitting matter densities 
f2m/i^ are consistent with the latest analyses of the CMB 
(Komatsu et al. 2011). The damping parameters cj^ are not 
well-constrained using our data, but the allowed range is 
consistent with the predictions of Equation 5 for our fidu- 
cial model (which are (j„ = (4.8, 4.5, 4.2) Mpc for the 
three redshift slices). When fitting (t„ we only permit it to 
vary over the range ai, > 0. 

The 2D probability contours for and a, marginal- 

izing over and are displayed in Figure 3. The measure- 
ment of a (hence Dy — a Dv,{id) is significantly correlated 
with the matter density, which controls the shape of the 
clustering pattern. 

We indicate three degeneracy directions in the param- 
eter space of Figure 3. The first direction (the dashed line) 
corresponds to a constant measured acoustic peak separa- 
tion, i.e. a/rs{zd) = constant, where rs{zd) is the sound 
horizon at the drag epoch as a function of Q,mh^ , deter- 
mined using the fitting formula quoted in Equation 12 of 
Percival et al. (2010). This parameter degeneracy would be 
expected in the case that just the baryon acoustic peak is 
driving the model fits, such that the measured low-redshift 
distance aDyfid is proportional to the standard ruler scale 
rs{zd)- 
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Figure 2. The amplitude of the cross-correlation Cij / ^CuCjj of the covariance matrix dj for the combined WiggleZ correlation 
function measurements for the redshift ranges 0.2 < 2 < 0.6, 0.4 < 2 < 0.8, 0.6 < z < 1.0 and 0.2 < 2 < 1.0, determined using lognormal 
realizations. 



The second direction (the dotted line) illustrated in Fig- 
ure 3 represents the degeneracy resulting from a constant 
measured shape of a Cold Dark Matter (CDM) power spec- 
trum, i.e. fluih^ X a = constant. We note here the consistency 
between this scaling and the "shape parameter" F = flmh 
used to parameterize the CDM transfer function (Bardeen et 
al. 1986). This shape parameter assumes that wavenumbers 
are observed in units of h Mpc~^, but the standard ruler 
scale encoded in baryon acoustic oscillations is calibrated 
by the CMB in units of Mpc, with no factor of h. 

The third direction (the dash-dotted line) shown in Fig- 
ure 3, which best describes the degeneracy in our data, cor- 
responds to a constant value of the acoustic parameter A(z) 
introduced by Eisenstein et al. (2005), 



A{z) 



100 Dviz)V^hJ^ 



(14) 



which appears in Figure 3 as \/jXn7? x q = constant. We 
note that the values of A{z) predicted by any cosmological 
model are independent of h, because Dv is proportional to 

h-\ 

The acoustic parameter A{z) provides the most ap- 
propriate description of the distance-redshift relation deter- 
mined by a BAG measurement in which both the clustering 
shape and acoustic peak are contributing toward the fit, 
such that the whole correlation function is being used as a 
standard ruler (Eisenstein et al. 2005, Sanchez et al. 2008, 



Shoji et al. 2009). In this case, the resulting measurement 
of A{z) is approximately uncorrelated with Q,mh^- We re- 
peated our BAG fit to the WiggleZ correlation functions in 
redshift slices using the parameter set {A, fJm/i'^, cr„, 6^). The 
marginalized values of A{z) we obtained are quoted in Table 
1, and correspond to measurements of the acoustic parame- 
ter with accuracies 7.2%, 4.5% and 5.0% in the three redshift 
slices. 

We also fitted our data with the parameter set 
(d^, fim/i^, cTu, 6^), where dz = rs{zd)/Dv{z). Results are 
again listed in Table 1, corresponding to measurements of 
dz with accuracies 7.8%, 4.7% and 5.4% in the three red- 
shift slices. We note that, unlike for the case of A{z), these 
measurements of dz are correlated with the matter density 
Quih^ , due to the orientation of the parameter degeneracy di- 
rections in Figure 3 (noting that constant dz corresponds to 
the "constant measured acoustic peak" case defined above). 

As a check for systematic modelling errors, we repeated 
the fits to the WiggleZ correlation functions using the second 
acoustic peak model described in Section 3, motivated by 
perturbation theory, fitting the data over the same range of 
scales. The marginalized measurements of a in the three red- 
shift slices were (1.032 ± 0.093, 0.981 ± 0.060, 1.091 ± 0.079), 
to be compared with the results for the default model quoted 
in Table 1. The amplitude of the systematic error in the fit- 
ted scale distortion parameter is hence significantly lower 
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Table 1. Results of fitting the four-parameter model (Hm/i^, a, cr„, fe'^) to the WiggleZ correlation functions in three redshift slices, 
together with the results for the full sample. The effective redshifts of the measurement in each slice are listed in Column 2, and the 
corresponding values of Dy for the fiducial cosmological model appear in Column 3. The values of for the best-fitting models are 
quoted in Column 4, for 13 degrees of freedom. Columns 5, 6 and 7 show the marginalized measurements of the matter density parameter 
flmh^, scale distortion parameter a and damping scale cr^ in each redshift slice. Corresponding measurements of the BAO distilled 
parameters A{z) and dz are displayed in Columns 8 and 9. The measured values of Dy in each redshift slice are given by aDv.fid- 



Sample 






Zcff 


Dv,iid 






Q 




A{z^b) 












[Mpc] 








Mpc] 






WiggleZ - 0.2 < 


z 


< 0.6 


0.44 


1617.8 


11.4 


0.143 ± 0.020 


1.024 ±0.093 


4.5 ±3.5 


0.474 ± 0.034 


0.0916 ±0.0071 


WiggleZ - 0.4 < 


z 


< 0.8 


0.60 


2085.4 


10.1 


0.147 ±0.016 


1.003 ±0.065 


4.1 ±3.4 


0.442 ± 0.020 


0.0726 ± 0.0034 


WiggleZ - 0.6 < 


z 


< 1.0 


0.73 


2421.9 


13.7 


0.120 ±0.013 


1.113 ±0.071 


4.4 ±3.2 


0.424 ± 0.021 


0.0592 ± 0.0032 


WiggleZ - 0.2 < 


z 


< 1.0 


0.60 


2085.4 


11.5 


0.127 ±0.011 


1.071 ±0.053 


4.4 ±3.3 


0.441 ±0.017 


0.0702 ± 0.0032 



than the statistical error in the measurement (by at least a 
factor of 3 in all cases). 

We assessed the statistical significance of the BAO de- 
tections in each redshift slice by repeating the parameter fits 
replacing the model correlation function with one generated 
using the "no- wiggles" reference power spectrum Prcf(fc) 
as a function of (Eisenstein & Hu 1998). The min- 

imum values obtained for the statistic for the fits in the 
three redshift slices were 15.2, 15.1 and 19.4, indicating that 
the model containing baryon oscillations was favoured by 
Ax^ = 3.8, 5.0 and 5.7 (with the same number of parame- 
ters fitted). These intervals correspond to detections of the 
baryon acoustic peaks in the redshift slices with statisti- 
cal significances between 1.9-cr and 2.4-(t. We note that the 
marginalized uncertainty in the scale distortion parameter 
for the no- wiggles model fit degrades by a factor of between 
two and three compared to the fit to the full model, demon- 
strating that the acoustic peak is very important for estab- 
lishing the distance constraints from our measurements. 

We used the same approach to determine the statistical 
significance of the BAO detection in the full WiggleZ redshift 
span 0.2 < z < 1.0, after combining the correlation function 
measurements in the redshift slices 0.2 < z < 0.6 and 0.6 < 
z < 1.0. In this case the model containing baryon oscillations 
was favoured by Ax^ = 7.7, corresponding to a statistical 
significance of 2.8-a for the detection of the baryon acoustic 
peak. 



4.3 Covariances between redshift slices 

We used the ensemble of lognormal realizations to quantify 
the covariance between the BAO measurements in the three 
overlapping WiggleZ redshift slices. For each of the 400 log- 
normal realizations in every WiggleZ region, we measured 
correlation functions for the redshift ranges Azi = 0.2 < 
z < 0.6, Az2 = 0.4 < 2 < 0.8 and Az3 = 0.6 < z < 1.0 and 
combined these correlation functions for the different regions 
using inverse-variance weighting. We then fitted the default 
clustering model described in Section 3 to each of the 400 
combined correlation functions for the three redshift slices. 

Figure 4 displays the correlations between the 400 
marginalized values of the scale-distortion parameter a for 
every pair of redshift slices. As expected, significant corre- 
lations are found in the values of a obtained from fits to 
the overlapping redshift ranges (Azi, Az2) and {Az2, Azg), 
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Figure 3. Probability contours of the physical matter density 
Qmh^ and scale distortion parameter a obtained by fitting to 
the WiggleZ survey combined correlation function in four red- 
shift ranges 0.2 < z < 0.6, OA < z < 0.8, 0.6 < z < 1.0 and 
0.2 < z < 1.0. The heavy dashed and dotted lines are the degen- 
eracy directions which are expected to result from fits involving 
respectively just the acoustic peak, and just the shape of a pure 
CDM power spectrum. The heavy dash-dotted line represents a 
constant value of the acoustic "A" parameter defined by Equa- 
tion 14, which is the parameter best-measured by the WiggleZ 
correlation function data. The solid circle represents the location 
of our fiducial cosmological model. The contour level in each case 
encloses regions containing 68.27% of the total likelihood. 

whereas the fits to the non-overlapping pair (AzijAzs) 
produce an uncorrelated measurement (within the statis- 
tical noise). The corresponding correlation coefficients for 
the overlapping pairs are pi2 = 0.369 and p23 = 0.438, 
where pij = Cij / y^CuCjj in terms of the covariances 
Cij = {aiOj} — {ai){aj). Table 2 contains the resulting in- 
verse covariance matrix for the measurements of A{z) in the 
three redshift slices, that should be used in cosmological pa- 
rameter fits. 

4.4 Comparison to mock galaxy catalogue 

As a further test for systematic errors in our distance scale 
measurements we fitted our BAO models to a dark matter 
halo catalogue generated as part of the Gigaparsec WiggleZ 
(GiggleZ) simulation suite (Poole et al., in prep.). The main 
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Figure 4. These panels illustrate the correlations between the scale distortion parameters a fitted to correlation functions for three 
overlapping WiggleZ redshift slices using 400 lognormal realizations. The red ellipses represent the derived correlation coefficients between 
these measurements. 



Table 2. The inverse covariance matrix C ^ of the measurements 

from the WiggleZ survey data of the acoustic parameter A{z) de- 
fined by Equation 14. We have performed these measurements in 
three overlapping redshift slices 0.2 < z < 0.6, 0.4 < z < 0.8 and 
0.6 < 2 < 1.0 with effective redshifts z^g = 0.44, 0.6 and 0.73, 
respectively. The data vector is A^^^, = (0.474,0.442,0.424), as 
listed in Table 1. The chi-squared statistic for any cosmological 
model vector A„^^ can be obtained via the matrix multiplication 



■Amod)'^Q. Miiobs -iimod)- The matrix is symmet- 



ric; we just quote the upper diagonal. 



Redshift slice 


0.2 < 2 < 0.6 


0.4 < 2 < 0.8 


0.6 < 2 < 1.0 


0.2 < 2 < 0.6 


1040.3 


-807.5 


336.8 


0.4 < 2 < 0.8 




3720.3 


-1551.9 


0.6 < 2 < 1.0 






2914.9 



GiggleZ simulation consists of a 2160 particle dark matter 
N-body calculation in a box of side 1 h^^ Gpc. The cosmo- 
logical parameters used for the simulation initial conditions 
were [n^, f^b, rie, /i, erg] = [0.273,0.0456,0.96,0.705,0.812]. 

We measured the redshift-space correlation function of 
a mass-limited subset of the dark matter halo catalogue ex- 
tracted from the z = 0.6 snapshot. This subset of dark mat- 
ter haloes, spanning a small range of maximum circular ve- 
locities around 125 km/s, was selected to possess a similar 
large-scale clustering amplitude to the WiggleZ galaxies at 
that redshift. We obtained the covariance matrix of the mea- 
surement using jack-knife techniques. We fitted our default 
correlation function model described in Section 3 to the re- 
sult, varying Qnih^ , a, a„ and and using the same fitting 
range as the WiggleZ measurement, 10 < s < 180 Mpc. 

Figure 5 shows the z — 0.6 GiggleZ halo correlation 
function measurement compared to the WiggleZ correlation 
function for the redshift range 0.4 < z < 0.8 (which was plot- 
ted in the top right-hand panel of Figure 1). We overplot the 
best-fitting default correlation function model for the Gig- 
gleZ data. The 2D probability contours for Qmh^ and a are 
displayed in Figure 6, again compared to the 0.4 < z < 0.8 
WiggleZ measurement and indicating the same degeneracy 
directions as shown in Figure 3. We conclude that the best- 
fitting parameter values are consistent with the input values 
of the simulation (within the statistical error expected in a 



{ WiggleZ 0.4<z<0.8 
^ GiggleZ halo catalogue 




50 100 150 
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Figure 5. Measurement of the galaxy correlation function l;{s) 
from a GiggleZ redshift-space halo subset at 2 = 0.6, chosen to 
possess a similar large-scale clustering amplitude to the WiggleZ 
galaxies at that redshift. We plot the correlation function in the 
combination ^(s) where s is the co-moving redshift-space sep- 
aration, and compare the result to the WiggleZ correlation func- 
tion for the redshift range 0.4 < 2 < 0.8. The best-fitting cluster- 
ing model to the GiggleZ measurement, varying the parameters 
Qmh^, ct, o-y and as described in Section 3, is overplotted as 
the solid line. 



measurement that uses a single realization) and there is no 
evidence for significant systematic error. We note that the ef- 
fective volume of the halo catalogue is slightly greater than 
that of the WiggleZ survey redshift range 0.4 < z < 0.8, 
hence the BAO measurements are more accurate in the case 
of GiggleZ. 



5 BARYON ACOUSTIC PEAK 

MEASUREMENT FROM THE FULL SLOAN 
DIGITAL SKY SURVEY LUMINOUS RED 
GALAXY SAMPLE 

In this Section we measure and fit the correlation function 
of the SDSS-LRG DR7-Full sample. This analysis is similar 
to that performed by Kazin et al. (2010a) for quasi-volume- 
limited sub-samples with z < 0.36, but now extended to a 
higher maximum redshift z = 0.44. We note that we assume 
a fiducial cosmology Q,n-i = 0.25 for this analysis, motivated 
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Figure 6. Probability contours of the physical matter density 
Q,uih? and scale distortion parameter a obtained by fitting the 
default correlation function model to the GiggleZ halo subset at 
2 = 0.6. We compare the result to the WiggleZ measurement in 
the redshift range 0.4 < z < 0.8 and overplot the same degener- 
acy directions as shown in Figure 3. The solid circle represents 
the location of our fiducial cosmological model. The contour level 
in each case encloses regions containing 68.27% of the total like- 
lihood. 

by the cosmological parameters used in the LasDamas sim- 
ulations (which we use to determine the covariance matrix 
of the measurement as described below in Section 5.2). The 
choice instead of Sim = 0.27, as used for the 6dFGS and 
WiggleZ analyses, would yield very similar results because 
the Alcock-Paczynski distortion between these cases is neg- 
ligible compared to the statistical errors in a. 

5.1 Correlation function measurement 

We measured the correlation function of the SDSS-LRG 
DR7-Full sample by applying the estimator of Equation 
11, using random catalogues constructed in the manner de- 
scribed in detail by Kazin et al. (2010a). For the purposes 
of the model fits in this Section we used separation bins of 
width 6.6 Mpc spanning the range 40 < s < 2QQh~^ 
Mpc, although we also determined results in IQh^^ Mpc 
bins in order to combine with the 6dFGS and WiggleZ cor- 
relation functions in Section 6 below. The measurement of 
the DR7-Full correlation function in 6.6 Mpc bins is 
displayed in the left-hand panel of Figure 7, where the er- 
ror bars are determined from the diagonal elements of the 
covariance matrix of 160 mock realizations, generated as de- 
scribed below in Section 5.2. The solid and dashed lines 
in Figure 7 are two best-fitting models, determined as ex- 
plained below in Section 5.3. 

The correlation function measurements in the separa- 
tion range 120 < s < 190 Mpc are higher than ex- 
pected in the best-fitting model. However, it is important 
to remember that these data points are correlated. The re- 
duced chi-squared statistics corresponding to these models 
are x^/dof = 1.1 — 1.2 (for 22 degrees of freedom), which 
fall well within the distribution of found in individual 
fits to the 160 mock catalogues, as shown in the right-hand 
inset in Figure 7. Kazin et al. (2010a) discussed the excess 
clustering measurement in SDSS-LRG subsamples and sug- 



gested that this is likely to result from sample variance. This 
is now reinforced by the fact that the independent-volume 
measurements from the WiggleZ and 6dFGS samples do not 
show similar trends of excess (see Figure 8). 

A potential cause of the stronger-than-expected clus- 
tering of LRGs on large scales is the effect of not mask- 
ing faint stars on random-catalogue generation. Ross et al. 
(2011) showed that apparent excess large-scale angular clus- 
tering measured in photometric LRG samples (Blake et al. 
2007, Padmanabhan et al. 2007, Thomas et al. 2011) is a sys- 
tematic effect imprinted by anti-correlations between faint 
stars and the galaxies, that can be corrected for by masking 
out regions around the stars. However, in the sparser SDSS- 
DR7 LRG sample the faint stars are uncorrelated with the 
galaxies at the angles of interest and do not introduce signif- 
icant systematic errors in the measured correlation function 
(A.Sanchez, private communication). 

5.2 LasDamas mock galaxy catalogues 

We simulated the SDSS-LRG correlation function measure- 
ment and determined its covariance matrix using the mock 
galax;y catalogues provided by the Large Suite of Dark Mat- 
ter Simulations (LasDamas, McBride et al. in prep.). These 
N-body simulations were generated using cosmological pa- 
rameters consistent with the WMAP 5-year fits to the CMB 
fiuctuations (Komatsu et al. 2009): [fim, fib, Ws, /i, o"g] = 
[0.25,0.04,1.0,0.7,0.8]. 

The LasDamas collaboration generated realistic LRG 
mock catalogues* by placing galaxies inside dark matter ha- 
los using a Halo Occupation Distribution (HOD; Berlind & 
Weinberg 2002). The HOD parameters were chosen to repro- 
duce the observed galaxy number density as well as the pro- 
jected two-point correlation function Wp{r-p) of the SDSS- 
LRG sample at separations 0.3 < rp < 30 Mpc. We 
used a suite of 160 LRG mock catalogues constructed from 
light cone samples with a mean number density n ~ 10~* h? 
Mpc~^. Each DR7-Full mock catalogue covers the redshift 
range 0.16 < z < 0.44 and reproduces the SDSS angular 
mask, corresponding to a total volume \.2h~'^ Gpc^. The 
mock catalogues were subsampled to match the observed 
redshift distribution of the LRGs. 

5.3 Correlation function modelling 

We extracted the scale of the baryon acoustic feature in 
the DR7-Full correlation function measurement by fitting 
for the scale distortion parameter a relative to a template 
correlation function ^fld using Equation 1, fitting over the 
separation range 40 < s < 200 Mpc. Together with the 
two correlation function models already described in Section 
3, the availability of the suite of LasDamas mock catalogues 
allows us to add a third template to use as ^gd: the mock- 
mean correlation function ^moan of all 160 realizations, which 
includes effects due to the non-linear growth of structure, 
redshift-space distortions, galaxy bias, light-coning and the 
observed 3D mask. 

The best-fitting model taking = Cmoan, marginal- 
izing over the correlation function amplitude, is displayed 

t http:/ /Iss. phy.vanderbilt.edu/lasdamas/ 
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Figure 7. The left-hand plot displays our correlation function measurement for the SDSS-LRG DR7-Full sample over the separation 
range 40 < s < 200 ft," 1 Mpc (where the error bars are the diagonal elements of the covariance matrix determined from 160 mock 
realizations). The solid line is the best-fitting model based on the mock-mean correlation function ^mcam and the dashed line is the 
best- fitting analytic perturbation-theory model correlation function ^pt based on Equation 8. The arrows point to the most likely peak 
position according to each model, where the longer arrow corresponds to the ^mean result Spcaii = 102.2 ± 2.8 Mpc. In the right-hand 
inset the reduced chi-squared statistic x^/dof = 1.1 (1.2) using ^moan (Jpt) for 22 degrees of freedom is compared with a histogram of 
the results fitting to the 160 individual mock realizations. The left-hand inset compares the measurement of dz=o.3l4 to the distribution 
found from the mocks; the offset of the measured result is due to the fact that the fiducial matter density Qyn = 0.25 used to generate 
the mocks is a little lower than the current best fits to cosmological data. The right-hand plot shows the distribution amongst the 160 
mocks of the difference in the chi-squared statistic between a model containing the baryon acoustic peak and a featureless model. The 
3.4-fT detection of the baryon acoustic feature that we find in DR7-Full (Ax^ = 11.9) falls well within the distribution of values found 
by applying a similar analysis to the mock catalogues. 



as the solid line in Figure 7, corresponding to a = 1.045. 
The statistic of the best fit is 24.2 (for 22 degrees of 
freedom). The most likely baryon acoustic peak position 
(determined using the method of Kazin et al. 2010a) is 
Spcait ~ 102.2 ± 2.8 Mpc (represented by the large arrow 
in Figure 7), where the quoted error is based on the sam- 
ple variance determined by performing the same analysis on 
all 160 mock catalogues. The corresponding measurement of 
the distilled BAO parameter is d^=o.3i4 = 0.1239 ± 0.0033. 
The distribution of measurements of for the 160 mocks 
is shown as the left-hand inset in Figure 7. We do not ex- 
pect the SDSS result (vertical lines) to coincide with unity, 
because of the difference between the true and fiducial cos- 
mological parameters. 

As a comparison, we also fitted to these data the two 
correlation function models described in Section 3, param- 
eterized by (dz, Jlm/i^, cr„, 6^). The marginalized measure- 
ments of dz for the two models were 0.1265 ± 0.0048 and 
0.1272 ±0.0050, consistent with our determination based on 
the mock-mean correlation function (which effectively uses 
fixed values of Qmh^ and a^). 

Our best-fitting analytic perturbation-theory model due 
to Crocce & Scoccimarro (2008) is displayed as the red 
dashed line in the left-hand panel of Figure 7. In this model 
we find that the best-fitting value of Sp^ak is correlated with 



cr„, although such changes produce offsets smaller than the 
l-(j statistical error in a (represented by the grey region 
around the short arrows in Figure 7). 



5.4 Significance of detection of the SDSS-LRG 
baryon acoustic feature 

We assessed the statistical significance of the detection of 
the baryon acoustic peak in the SDSS-LRG sample in a 
similar style to the WiggleZ analysis described in Section 
4.2, by comparing the best-fitting values of x^ for models 
containing a baryon acoustic feature (Xfoaturc) ^nd feature- 
less models (Xfeatureiess) cousttucted using the "no- wiggles" 
power spectrum of Eisenstein & Hu (1998). We used the 
perturbation-theory model for the baryon acoustic peak de- 
scribed in Section 3 when constructing these models. 

The SDSS-LRG dataset produced Ax^ = xLturo - 
Xfcaturoicss ~ —11-9 ovcr the separation range 40 < s < 
200 Mpc, corresponding to a detection of the baryon 
acoustic feature with significance of SA-a. The histogram 
resulting from repeating this analysis for all 160 mocks is dis- 
played in the right-hand panel of Figure 7, following Cabre 
& Gaztanaga (2011); we see that the SDSS result is as ex- 
pected from an average realization. 

We used the same method to compare the signifi- 
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Figure 8. The correlation function measurements ^(s) for the WiggleZ, SDSS-LRG and 6dFGS galaxy samples, plotted in the combination 
(^(s) where s is the co- moving redshift-space separation. The lower right-hand panel shows the combination of these measuremements 
with inverse- variance weighting. The best-fitting clustering models in each case, varying the parameters f2ni^^, ct, ctv and 6^ as described 
in Section 3, are overplotted as the solid lines. 



cance of detection of the acoustic peak in DR7-Full with 
that obtained in the volume-limited LRG sub-samples an- 
alyzed by Kazin et al. (2010a). The sample "DR7-Sub", 
a quasi-volume-limited LRG catalogue spanning redshift 
range 0.16 < 2: < 0.36 and luminosity range —23.2 < Mg < 
—21.2, yields a detection significance of 2.2-a. For the sam- 
ple "DR7-Bright" , a sparser volume-limited catalogue with 
a brighter luminosity cut —23.2 < Mg < —21.8, the signifi- 
cance of the baryon acoustic feature is just below 2-a. 



6 THE STACKED BARYON ACOUSTIC PEAK 

Our goal in this Section is to assess the overall statistical 
significance with which the baryon acoustic peak is detected 
in the combination of current galaxy surveys. In order to do 
this we combined the galaxy correlation functions measured 
from the WiggleZ Survey, the Sloan Digital Sky Survey Lu- 
minous Red Galaxy (SDSS-LRG) sample and the 6-degree 
Field Galaxy Survey (6dFGS), and fitted the models de- 
scribed in Section 3 to the result. Although we acknowledge 
that model fits to a combination of correlation functions ob- 
tained using different redshifts and galaxy types will produce 
parameter values that evade an easy physical interpretation, 
the resulting statistical significance of the BAO detection re- 
mains a quantity of interest. 



6.1 The 6dFGS baryon acoustic peak 
measurement 

For completeness we summarize here the measurement of the 
baryon acoustic peak from the 6dFGS reported by Rentier 
et al. (2011). After optimal weighting of the data to min- 
imize the correlation function error at the baryon acous- 
tic peak, the 6dFGS sample covered an effective volume 
Kff = 0.08 /i"^ Gpc^ with effective redshift z^a = 0.106. 
Rentier et al. fitted the model defined by our Equation 8 
to the 6dFGS correlation function, using lognormal real- 
izations to determine the data covariance matrix and vary- 
ing the parameter set Q,mh^, a, o-y and b'^ . The model fits 
were performed over the separation range 10 < s < 190 
Mpc, with checks made that the best-fitting parameters were 
not sensitive to the minimum separation employed. The re- 
sulting measurements of the distance scale were quantified 
as DviO.lOG) = 457 ± 27 Mpc, do.ioe = 0.336 ± 0.015 or 
A(0.106) = 0.526 ± 0.028. The statistical significance of the 
detection of the acoustic peak was estimated to be 2.4-a, 
based on the difference in chi-squared Ax^ = 5.6 between 
the best-fitting model and the corresponding best fit of a 
zero-baryon model. 

6.2 The combined correlation function 

Figure 8 displays the three survey correlation functions 
combined in our study: the WiggleZ 0.2 < z < 1.0 mea- 
surement plotted in the lower right-hand panel of Figure 
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Figure 9. The amplitude of the cross-correlation dj / -^J Cu Cjj of the covariance matrix dj for the WiggleZ, SDSS-LRG and 6dFGS 
correlation functions. The lower right-hand panel shows the covariance matrix of the combined correlation function. The covariance 
matrices for the WiggleZ and 6dFGS samples are determined using lognormal realizations, and that of the SDSS-LRG sample is obtained 
from an ensemble of N-body simulations. The plot of the WiggleZ cross-correlation matrix in the upper-left hand panel is reproduced 
from the lower right-hand panel of Figure 2. 



1, the 6dFGS correlation function reported by Beutler et 
al. (2011), and the SDSS-LRG DR7-Full measurement de- 
scribed in Section 5 (using a binning of IQh'^ Mpc in aU 
cases). These correlation functions have quite different am- 
plitudes owing to differences between the growth factors at 
the effective redshifts z of the samples and the bias factors 
b of the various galaxy tracers. Before stacking these func- 
tions we make an amplitude correction to a common red- 
shift zo = 0.35 and bias factor bo — 1, by multiplying each 
correlation function by [bl G{zof Bo{l3o)]/[b^ G{zf Bo{l3)] 
where G{z) is the linear growth factor at redshift z and 
Bo{P) — l-f|/3-|-i/3^ is the Kaiser boost factor in terms 



of the redshift-space distortion parameter /3 = 



N6/11 



/b 



(Kaiser 1987). When calculating these quantities we as- 
sumed that the redshifts of the WiggleZ, SDSS-LRG and 
6dFGS samples were z = (0.6, 0.314, 0.106) and the bias fac- 
tors were b = (1.1,2.2, 1.8). After making these normaliza- 
tion corrections we then combined the correlation functions 
and their corresponding covariance matrices using inverse- 
variance weighting in the same style as Equations 12 and 
13. The resulting total correlation function is plotted in the 
lower right-hand panel of Figure 8. The covariance matrices 
of the different survey correlation functions and final combi- 
nation are displayed in Figure 9. An additional overplot of 
these correlation functions is provided in Figure 10. We note 



that although the SDSS-LRG correlation function measure- 
ment used the fiducial cosmology Qm ~ 0.25, compared to 
the choice fim ~ 0.27 for the WiggleZ and 6dFGS analy- 
ses, the Alcock-Paczynski distortion between these cases is 
negligible compared to the statistical errors in a. 



6.3 Significance of the detection of the baryon 
acoustic peak in the combined sample 

We fitted the clustering model described in Section 3 to 
the combined correlation function over separation range 
30 < s < 180 Mpc, varying Qmh^, a, and b^ and 
using an effective redshift z = 0.35. We used the more con- 
servative minimum fitted scale 30 Mpc for the analysis of 
the stacked correlation function in this Section, compared to 
10 Mpc for the fits to the WiggleZ correlation function in 
Section 4, because (1) the required non- linear corrections be- 
come more important for galaxy samples such as the 6dFGS 
and SDSS LRGs, which are both more biased and at lower 
redshift than the WiggleZ sample, and (2) systematic errors 
in the fitting become relatively more important for this com- 
bined dataset with higher signal-to-noise. Although we fixed 
the relative bias factors of the galaxy samples when stack- 
ing the survey correlation functions in Section 6.2, we still 
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Figure 10. An overplot of the correlation function measure- 
ments e,{s) for the WiggleZ, SDSS-LRG and 6dFGS galaxy sam- 
ples, plotted in the combination ^(s) where s is the co-moving 
redshift-space separation. A normalization correction has been 
applied to these correlation functions to allow for the differing 
effective redshifts and galaxy bias factors of the samples (see text 
for details). The combined correlation function, determined by 
inverse-variance weighting, is also plotted. The best-fitting clus- 
tering model to the combined correlation function (varying Q.uih? , 
o, cr„ and b^) is overplotted as the solid line. We also show as 
the dashed line the corresponding "no-wiggles" reference model 
(Eisenstein &; Hu 1998), constructed from a power spectrum with 
the same clustering amplitude but lacking baryon acoustic oscil- 
lations. 



marginalized over an absolute normalization 6^ ~ 1 when 
fitting the model in this Section. 

We obtained a good fit to the stacked correlation func- 
tion with = 11-3 (for 11 degrees of freedom) and 
marginalized parameter values Q.^h? = 0.132 ± 0.014, a = 
1.037±0.036 ando-t, = 4.5±1.8/i"^ Mpc. Although the best- 
fitting value of a must be interpreted as some effective value 
integrating over redshift, we can conclude that the measured 
BAO distance scale is consistent with the fiducial model. 

We quantified the significance of the detection of the 
acoustic peak in the combined sample using two methods. 
Firstly, we repeated the parameter fit replacing the model 
correlation function with one generated using a "no-wiggles" 
reference power spectrum (Eisenstein & Hu 1998). The mini- 
mum value obtained for the statistic in this case was 32.7, 
indicating that the model containing baryon oscillations was 
favoured by A^^ = 21.4. This corresponds to a detection of 
the acoustic peak with a statistical significance of 4.6-(t. 

As an alternative approach for assessing the significance 
of the detection, we changed the fiducial baryon density to 
f2b = and repeated the parameter fit. For zero baryon den- 
sity we generated the model matter power spectrum using 
the fitting formulae of Eisenstein & Hu (1998), rather than 
using the CAMB software. The minimum value obtained for 
the statistic was now 35.3, this time suggesting that the 
acoustic peak had been detected with a significance of 4.9-cr. 
The reason that the significance of detection varies between 
these two methods of assessment is that in the latter case, 
where the baryon density is changed, the overall shape of 
the clustering pattern is also providing information used to 
disfavour the S7b = model, whereas in the former case only 



the presence of the acoustic peak varies between the two sets 
of models. 



7 COSMOLOGICAL PARAMETER FITS 

In this Section we fit cosmological models to the latest 
distance datasets comprising BAO, supernovae and CMB 
measurements. Our aim is to compare parameter fits to 
BAO-f-CMB data (excluding supernovae) and SNe-|-CMB 
data (excluding BAO) as a robust check for systematic er- 
rors in these distance probes. 



7.1 BAO dataset 

The latest BAO distance dataset, including the 6dFGS, 
SDSS and WiggleZ surveys, now comprises BAO mea- 
surements at six different redshifts. These data are sum- 
marized in Table 3. Firstly, we use the measurement of 
do. 106 = 0.336 ± 0.015 from the 6dFGS reported by Beutler 
et al. (2011). Secondly, we add the two correlated measure- 
ments of do. 2 and do. 35 determined by Percival et al. (2010) 
from fits to the power spectra of LRGs and main-sample 
galaxies in the SDSS (spanning a range of wavenumbers 
0.02 < fc < 0.3/1 Mpc~^). The correlation coefficient for 
these last two measurements is 0.337. We note that our own 
LRG baryon acoustic peak measurements reported above in 
Section 5 are entirely consistent with these fits. Finally, we 
include the three correlated measurements of A{z = 0.44), 
A{z = 0.6) and A{z — 0.73) reported in this study, using 
the inverse covariance matrix listed in Table 2. 

In our cosmological model fitting we assume that the 
BAO distance errors are Gaussian in nature. Modelling po- 
tential non-Gaussian tails in the likelihood is beyond the 
scope of this paper, although we note that they may not 
be negligible (Percival et al. 2007, Percival et al. 2010, Bas- 
sett & Afshordi 2010). We caution that the 2-a confidence 
regions displayed in the Figures in this Section might not 
necessarily follow the Gaussian scaling. The WiggleZ and 
SDSS-LRG surveys share a sky overlap of ~ 500 deg^ for 
redshift range z < 0.5; given that the SDSS-LRG measure- 
ment is derived across a sky area ~ 8000 deg'^ and the errors 
in both measurements contain a significant component due 
to shot noise, the resulting covariance is negligible. 

This BAO distance dataset is plotted in Figure 11 rela- 
tive to a flat ACDM cosmological model with matter density 
fim ~ 0.29 and Hubble parameter h = 0.69 (these values 
provide the best fit to the combined cosmological datasets 
as discussed below). The panels of Figure 11 show various 
representations of the BAO dataset including Dy (z) and the 
distilled parameters A{z) and dz- 



7.2 SNe dataset 

We used the "Union 2" compilation by AmanuUah et al. 
(2010) as our supernova dataset, obtained from the website 
http://supernova.lbl.gov/Union. This compilation of 557 
supernovae includes data from Hamuy et al. (1996), Riess 
et al. (1999, 2007), Astier et al. (2006), Jha et al. (2006), 
Wood-Vasey et al. (2007), Holtzman et al. (2008), Hicken et 
al. (2009) and Kessler et al. (2009). The data is represented 
as a set of values of the distance modulus for each supernova 
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Figure 11. Current measurements of the cosmic distance scale using the BAO standard ruler applied to the 6dFGS, SDSS and WiggleZ 
surveys (where the data is taken from Beutler et al. 2011, Percival et al. 2010 and this study). The results are compared to a flat ACDM 
cosmological model with matter density f2m = 0.29 and Hubble parameter h = 0.69. Various representations of the data are shown: the 
BAO distance Dy{z) recovered from fits to the angle-averaged clustering measurements (top left-hand panel), these distances ratioed to 
the fiducial model (top right-hand panel), the distilled parameter A{z) (defined by Equation 14) extracted from fits governed by both 
the acoustic peak and clustering shape (bottom left-hand panel), and the distilled parameter dz determined by fits controlled by solely 
the acoustic peak information (bottom right-hand panel). We note that the conversion of the BAO fits to the measurements of Dv{z) 
presented in the upper two plots requires a value for the standard ruler scale to be assumed: we take ra(zrf) = 152.40 Mpc, obtained 
using Equation 6 in Eisenstein & Hu (1998) evaluated for our fiducial model Umh' = 0.1381 and Q^h'^ = 0.02227. 



Table 3. The BAO distance dataset from the 6dFGS, SDSS and 
WiggleZ surveys. Measurements of the distilled parameters dz 
and A(z) are quoted. The most appropriate choices to be used 
in cosmological parameter fits are indicated by bold font. For the 
SDSS data, the values of A{z) are obtained by scaling from the 
measurements of d^ reported by Percival et al. (2010) using their 
fiducial cosmological parameters and the same fractional error. 
The pairs of measurements at 2 = (0.2,0.35), z = (0.44,0.6) 
and z = (0.6, 0.73) are correlated with coefficients 0.337, 0.369 
and 0.438, respectively. The inverse covariance matrix of the data 
points at z = (0.2, 0.35) is given by Equation 5 in Percival et al. 
(2010). The inverse covariance matrix of the data points at z = 
(0.44, 0.6, 0.73) is given in Table 2 above. The other measurements 
are uncorrelated. 
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z 




A{z) 
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0.526 ± 0.028 


SDSS 


0.2 


0.1905 ±0.0061 
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SDSS 
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0.1097 ±0.0036 
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WiggleZ 
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WiggleZ 


0.6 


0.0726 ± 0.0034 


0.442 ±0.020 


WiggleZ 


0.73 


0.0592 ± 0.0032 


0.424 ± 0.021 



IMpc 



+ 25, 
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where Dl{z) is the luminosity distance at redshift z. The 
values of ^ are reported for a particular choice of the normal- 
ization Af — 51ogjQ/i, which is marginalized as an unknown 
parameter in our analysis as described below. When fitting 
cosmological models to these SNe data we used the full co- 
variance matrix of these measurements including systematic 
errors, as reported by AmanuUah et al. (2010). 

Figure 12 is a representation of the consistency and rel- 
ative accuracy with which baryon oscillation measurements 
and supernovae currently map out the cosmic distance scale. 
In order to construct this figure we converted the BAO mea- 
surements of Dv{z) into Da{z) assuming a Hubble param- 
eter H{z) for a flat ACDM model with flm = 0.29 and 
h — 0.69. The binned supernovae data currently measure 
the distance-redshift relation at z < 0.8 with 3 — 4 times 
higher accuracy than the BAOs, although we note that the 
consequences for cosmological parameter fits are highly in- 
fluenced by the differing normalization of the two methods. 
The supernovae measure the relative luminosity distance to 
the relation at z = 0, Dl{z)Ho/c, owing to the unknown 
value of the standard-candle absolute magnitude M. The 
BAOs measure a distance scale relative to the sound hori- 
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Figure 12. Comparison of the accuracy with which supernovae 
and baryon acoustic oscillations map out the cosmic distance 
scale at z < 0.8. For the purposes of this Figure, BAO mea- 
surements of Dv{z) have been converted into Da{z) assum- 
ing a Hubble parameter H{z) for a flat ACDM model with 
f2in = 0.29 and h = 0.69, indicated by the solid line in the Fig- 
ure, and SNe measurements of Di^(z) have been plotted assuming 
Da{z) = Dl(z)/{1 + zY. 



zon at baryon drag calibrated by the CMB data, effectively 
an absolute measurement of Dv{z) given that the error is 
dominated by the statistical uncertainty in the clustering 
fits, rather than any systematic uncertainty in the sound 
horizon calibration from the CMB. 

When undertaking cosmological fits to the supernovae 
dataset, we performed an analytic marginalization over the 
unknown absolute normalization A/ — 51ogj^Q/i (Goliath et al. 
2001, Bridle et al. 2002). This is carried out by determining 
the chi-squared statistic for each cosmological model as 



2 T ^-1 

X Q. y 

~ SN ~ 



(16) 



where y is the vector representing the difference between the 
distance moduli of the data and model, and C_~^ is the in- 

SN 

verse covariance matrix for the supernovae distance moduli. 



7.3 CMB dataset 

We included the CMB data in our cosmological fits using 
the Wilkinson Microwave Anisotropy Probe (WMAP) "dis- 
tance priors" (Komatsu et al. 2009) using the 7-year WMAP 
results reported by Komatsu et al. (2011). The distance pri- 
ors quantify the complete CMB likelihood via a 3-parameter 
covariance matrix for the acoustic index Ia, the shift param- 
eter TZ and the redshift of recombination z,, as given in Table 
10 of Komatsu et al. (2011). When deriving these quantities 
we assumed a physical baryon density fly^h^ = 0.02227, a 
CMB temperature Tcmb ~ 2.72^K and a number of rela- 
tivistic degrees of freedom A'cft = 3.04. 



7.4 Flat vj models 

We first fitted a flat luCDM cosmological model in which 
spatial curvature is fixed at f^k = but the equation-of-state 



Figure 13. The joint probability for parameters S^m and w fitted 
separately to the WMAP, BAO and SNe distance data, marginal- 
ized over nni/i^ and assuming £7^ = 0. The two contour levels in 
each case enclose regions containing 68.27% and 95.45% of the 
total likelihood. 



w of dark energy is varied as a free parameter. We fitted 
for the three parameters (fim, fim/i^, w) using flat, wide pri- 
ors which extend well beyond the regions of high likelihood 
and have no effect on the cosmological fits. The best-fitting 
model has = 532.9 for 563 degrees of freedom, represent- 
ing a good fit to the distance dataset. 

Figures 13 and 14 compare the joint probability of fim 
and w, marginalizing over f2mft'^, for the individual WMAP, 
BAO and SNe datasets along with various combinations. 
We note that for the "BAO only" contours in Figure 13, we 
have not used any CMB calibration of the standard ruler 
scale rs{zd), and thus the 6dFGS and SDSS measurements of 
dz = rs{zd) /Dv{z) do not contribute strongly to these con- 
straints. Hence the addition of the CMB data in Figure 14 
has the benefit of both improving the information from the 
dz measurements by determining rs{zd), and contributing 
the WMAP distance prior constraints. The WMAP-I-BAO 
and WMAP-l-SNe data produce consistent determinations of 
the cosmological parameters, with the error in the equation- 
of-state Aw « 0.1. Combining all three datasets produces 
the marginalized result w = —1.034 ± 0.080 (errors in the 
other parameters are listed in Table 4; the quoted error in h 
results from fitting the three parameters fim, h and w). The 
best-fitting equation-of-state is consistent with a cosmolog- 
ical constant model for which w = —1. 

We caution that the probability contours plotted in Fig- 
ures 13 and 14 (and other similar Figures in this Section) as- 
sume that the errors in the BAO distance dataset are Gaus- 
sian. If the likelihood contains a significant non-Gaussian 
tail, the 2-a region could be affected. 

We repeated the WMAP-I-BAO fit comparing the two 
different implementations of the SDSS-LRG BAO distance- 
scale measurements: the Percival et al. (2010) power spec- 
trum fitting at z — 0.2 and z = 0.35, and our corre- 
lation function fit presented in Section 5. We found that 
the marginalized measurements of w in the two cases were 
— 1.00 ±0.13 and —0.97 ±0.13, respectively. Our results are 
therefore not significantly changed by the methodology used 
for these LRG fits. 
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Figure 14. The joint probability for parameters f2ni and w fitted 
to various combinations of WMAP, BAO and SNe distance data, 
marginalized over Qmh'^ and assuming Q,^ = 0. The two contour 
levels in each case enclose regions containing 68.27% and 95.45% 
of the total likelihood. 
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Figure 15. The joint probability for parameters f!m and fij^ fitted 
separately to the WMAP, BAO and SNe distance data, marginal- 
ized over Qm'i^ and assuming to = — f. The two contour levels 
in each case enclose regions containing 68.27% and 95.45% of the 
total likelihood. 



7.5 Curved A models 

We next fitted a curved ACDM model, in which we fixed the 
equation-of-state of dark energy at = — 1 but added the 
spatial curvature as an additional free parameter. We 
fitted for the three parameters (nm, f2m/i'^, Hk) using flat, 
wide priors which extend well beyond the regions of high 
likelihood and have no effect on the cosmological fits. The 
best-fitting model has = 532.7 for 563 degrees of freedom. 

Figures 15 and 16 compare the joint probability of 
fim and r^k, marginalizing over S7m^'^, for the individual 
WMAP, BAO and SNe datasets along with various com- 
binations. Once more, we find that fits to WMAP-I-BAO 
and WMAP -(-SNe produce mutually consistent results. The 
BAO data has higher sensitivity to curvature because of 
the long lever arm represented by the relation of distance 
measurements at z < 1 and at recombination. Combining 
all three datasets produces the marginalized result Hk = 
—0.0040 ± 0.0062 (errors in the other parameters are listed 
in Table 4). The best-fitting parameters are consistent with 
zero spatial curvature. 

7.6 Additional degrees of freedom 

We fitted two further cosmological models, each containing 
an additional parameter. Firstly we fitted a curved kjCDM 
model in which we varied both the dark energy equation- 
of-state and the spatial curvature as free parameters. The 
best-fitting model has = 531.9 for 562 degrees of free- 
dom, representing an improvement of Ax^ = 1-0 compared 
to the case where fik = 0, for the addition of a single ex- 
tra parameter. In terms of information criteria this does 
not represent a sufficient improvement to justify the addi- 
tion of the extra degree of freedom. Figure 17 compares the 
joint probability of lo and f^k, marginalizing over fim and 
for the three cases WMAP+BAO, WMAP+SNe and 
WMAP-l-BAO-l-SNe. Combining all three datasets produces 
the marginalized measurements w — —1.063 ± 0.094 and 
fik = -0.0061 ± 0.0070. 

We finally fitted a flat w(a)CDM cosmological model in 
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Figure 16. The joint probability for parameters Qm and Qk fitted 
to various combinations of WMAP, BAO and SNe distance data, 
marginalized over Q^h^ and assuming w = —1. The two contour 
levels in each case enclose regions containing 68.27% and 95.45% 
of the total likelihood. 

which spatial curvature is fixed at f2k = but the equation- 
of-state of dark energy is allowed to vary with scale factor 
a in accordance with the Chevallier-Polarski-Linder param- 
eterization 'w{a) = Wo + (1 — ci)wa (Chevallier & Polarski 
2001, Linder 2003). The best-fitting model has x^ = 531.9 
for 562 degrees of freedom, and again the improvement in the 
value of x^ compared to the case where Wa = does not jus- 
tify the addition of the extra degree of freedom. Combining 
all three datasets produces the marginalized measurements 
wo = -1.09 ± 0.17 and Wa = 0.19 ± 0.69. We note that 
the addition of the BAO measurements to the WMAP-|-SNe 
dataset produces a more significant improvement for fits in- 
volving f2k than for Wa- 

In all cases, the best-fitting parameters are consistent 
with a flat cosmological constant model for which wo = — 1, 
Wa = and Slk ~ 0. The best-fltting values and errors in 
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Figure 17. The joint probability for parameters and w fitted 
to various combinations of WMAP, BAO and SNe distance data, 
marginalized over f2m and Vl^h^ . The two contour levels in each 
case enclose regions containing 68.27% and 95.45% of the total 
likelihood. 
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Figure 18. The joint probability for parameters uiq and Wa de- 
scribing an evolving equation-of-state for dark energy, fitted to 
various combinations of WMAP, BAO and SNe distance data, 
marginalized over f!m and f7ni/i^ and assuming fijj = 0. The two 
contour levels in each case enclose regions containing 68.27% and 
95.45% of the total likelihood. 



the parameters for the various models, for the fits using all 
three datasets, are listed in Table 4. 



8 CONCLUSIONS 

We summarize the results of our study as follows: 

• The final dataset of the WiggleZ Dark Energy Survey 
allows the imprint of the baryon acoustic peak to be detected 
in the galaxy correlation function for independent redshift 
slices of width Az = 0.4. A simple quasi-linear acoustic peak 
model provides a good fit to the correlation functions over 
a range of separations 10 < s < 180 Mpc. The result- 
ing distance-scale measurements are determined by both the 
acoustic peak position and the overall shape of the clus- 
tering pattern, such that the whole correlation function is 



being used as a standard ruler. As such, the acoustic param- 
eter A[z) introduced by Eisenstein et al. (2005) represents 
the most appropriate distilled parameter for quantifying the 
WiggleZ BAO measurements, and we present in Table 2 a 
3x3 covariance matrix describing the determination of A{z) 
from WiggleZ data at the three redshifts z = 0.44, 0.6 and 
0.73. We test for systematics in this measurement by vary- 
ing the fitting range and implementation of the quasi-linear 
model, and also by repeating our fits for a dark matter halo 
subset of the Gigaparsec WiggleZ simulation. In no case do 
we find evidence for significant systematic error. 

• We present a new measurement of the baryon acoustic 
feature in the correlation function of the Sloan Digital Sky 
Survey Luminous Red Galaxy (SDSS-LRG) sample, finding 
that the feature is detected within a subset spanning the 
redshift range 0.16 < z < 0.44 with a statistical significance 
of 3.4-(T. We derive a measurement of the distilled parameter 
rfz=o.3i4 ~ 0.1239 ± 0.0033 that is consistent with previous 
analyses of the LRG power spectrum. 

• We combine the galaxy correlation functions measured 
from the WiggleZ, 6-degree Field Galaxy Survey and SDSS- 
LRG samples. Each of these datasets shows independent ev- 
idence for the baryon acoustic peak, and the combined cor- 
relation function contains a BAO detection with a statistical 
significance of 4.9-(t relative to a zero-baryon model with no 
peak. 

• We fit cosmological models to the combined 6dFGS, 
SDSS and WiggleZ BAO dataset, now comprising six 
distance-redshift data points, and compare the results to 
similar fits to the latest compilation of supernovae (SNe) 
and Cosmic Microwave Background (CMB) data. The BAO 
and SNe datasets produce consistent measurements of the 
equation-of-state w of dark energy, when separately com- 
bined with the CMB, providing a powerful check for sys- 
tematic errors in either of these distance probes. Combining 
all datasets, we determine w = — 1.034±0.080 for a flat Uni- 
verse, and ilk = — 0.0040±0.0062 for a curved, cosmological- 
constant Universe. 

• Adding extra degrees of freedom always produces best- 
fitting parameters consistent with a cosmological constant 
dark-energy model within a spatially-flat Universe. Vary- 
ing both curvature and w, we find marginalized errors 
w = -1.063 ± 0.094 and fik = -0.0061 ± 0.0070. For a 
dark-energy model evolving with scale factor a such that 
iu(a) = «;o -I- (1 — a)wa, we find that wo ~ —1.09 ± 0.17 and 
Wa =0.19 ±0.69. 



In conclusion, we have presented and analyzed the most 
comprehensive baryon acoustic oscillation dataset assembled 
to date. Results from the WiggleZ Dark Energy Survey have 
allowed us to extend this dataset up to redshift z = 0.73, 
thereby spanning the whole redshift range for which dark en- 
ergy is hypothesized to govern the cosmic expansion history. 
By fitting cosmological models to this dataset we have es- 
tablished that a flat ACDM cosmological model continues to 
provide a good and self-consistent description of CMB, BAO 
and SNe data. In particular, the BAO and SNe yield con- 
sistent measurements of the distance-redshift relation across 
the common redshift interval probed. Our results serve as a 
baseline for the analysis of future CMB datasets provided 
by the Planck satellite (Ade et al. 2011) and BAO mea- 
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Table 4. The results of fitting various cosmological models to a combination of the latest CMB, BAO and SNc distance datascts. 
Measurements and l-a errors are listed for each parameter, marginalizing over the other parameters of the model. All models contain 
either (Clm, ^mh'^) or (Clin,h) amongst the parameters fitted. 



Model 




d.o.f. 


r2ni 




h 




Wo 


Wa 


Flat ACDM 


533.1 


564 


0.290 ±0.014 


0.1382 ±0.0029 


0.690 ± 0.009 








Flat juCDM 


532.9 


563 


0.289 ±0.015 


0.1395 ±0.0043 


0.696 ± 0.017 




-1.034 ± 0.080 




Curved ACDM 


532.7 


563 


0.292 ±0.014 


0.1354 ±0.0054 


0.681 ±0.017 


-0.0040 ± 0.0062 






Curved toCDM 


531.9 


562 


0.289 ±0.015 


0.1361 ±0.0055 


0.687 ± 0.019 


-0.0061 ± 0.0070 


-1.063 ± 0.094 




Flat «;(a)CDM 


531.9 


562 


0.288 ±0.016 


0.1386 ±0.0053 


0.695 ± 0.017 




-1.094 ±0.171 


0.194 ±0.687 



surements from the Baryon Oscillation Spectroscopic Survey 
(Eisenstein et al. 2011). 



ACKNOWLEDGMENTS 

We thank the anonymous referee for careful and constructive 
comments that improved this study. 

Wc acknowledge financial support from the Aus- 
tralian Research Council through Discovery Project grants 
DP0772084 and DP1093738 funding the positions of SB, DP, 
MP, GP and TMD. SC and DC acknowledge the support of 
the Australian Research Council through QEII Fellowships. 
MJD thanks the Gregg Thompson Dark Energy Travel Fund 
for financial support. 

We thank the LasDamas project for making their mock 
catalogues publicly available. In particular EK is much 
obliged to Cameron McBride for supplying mock catalogues 
on demand. EK also thanks Ariel Sanchez for fruitful lengthy 
discussions. EK was partially supported by a Google Re- 
search Award and NASA Award. 

FB is supported by the Australian Government through 
the International Postgraduate Research Scholarship (IPRS) 
and by scholarships from ICRAR and the A AO. 

GALEX (the Galaxy Evolution Explorer) is a NASA 
Small Explorer, launched in April 2003. Wc gratefully ac- 
knowledge NASA's support for construction, operation and 
science analysis for the GALEX mission, developed in co- 
operation with the Centre National d'Etudes Spatiales of 
Prance and the Korean Ministry of Science and Technology. 

Finally, the WiggleZ survey would not be possible with- 
out the dedicated work of the staff of the Australian Astro- 
nomical Observatory in the development and support of the 
AAOmega spectrograph, and the running of the AAT. 



REFERENCES 

Ade P., et al., 2011, A&A accepted (arXiv: 1101. 2022) 
AmanuUah R., et al., 2010, ApJ, 716, 712 
Astier P., et al., 2006, A&A, 447, 31 

Bardeen J.M., Bond J.R., Kaiser N., Szalay A.S., 1986, ApJ, 304, 
15 

Bassett B., Afsliordi N., 2010 (arXiv : 1005 . 1664) 
Berlind A.A., Weinberg D.H., 2002, ApJ, 575, 587 
Beutlcr F., et al., 2011, MNRAS accepted (arXiv: 1106.3366) 
Blake C.A, Glazebrook K., 2003, ApJ, 594, 665 
Blake C.A., Parkinson D., Bassett B., Glazebrook K., Kunz M., 
Nichol R.C., 2006, MNRAS, 365, 255 



Blake C.A., CoUister A., Bridle S., Lahav O., 2007, MNRAS, 374, 
1527 

Blake C.A., et al., 2011, MNRAS accepted (arXiv: 1105. 2862) 

Bond J.R., Efstathiou G., 1984, ApJ, 285, L45 

Bridle S.L., Crittenden R., Melchiorri A., Hobson M.P., Kneissl 

R., Lascnby A.N., 2002, MNRAS, 335, 1193 
Cabrc A., Gaztanaga E., 2011, MNRAS, 412, 98 
Chevallier M., Polarski D., 2001, Int. J. Mod. Phys., DIO, 213 
Cole S., et al., 2005, MNRAS, 362, 505 

Cooray A., Hu W., Huterer D., Joffre M., 2001, ApJ, 557L, 7 
Copin Y., et al., 2006, New Astronomy Review, 50, 436 
Crocce M., Scoccimarro R., 2008, PRD, 2008, 77, 3533 
Crocce M., Gaztanaga E., Cabre A., Carnero E., Sanchez E., 2011, 

MNRAS accepted (arXiv: 1104.5236) 
Dawson K.S., et al., 2009, AJ, 138, 1271 
Drinkwater M., et al., 2010, MNRAS, 401, 1429 
Efstathiou G., Sutherland W., Maddox S., 1990, Nature, 348, 705 
Eisenstein D.J., Hu W., 1998, ApJ, 496, 605 
Eisenstein D.J., Hu W., Tegmark M., 1998, ApJ, 504, 57 
Eisenstein D.J., et al., 2001, AJ, 122, 2267 

Eisenstein D.J., 2003, "Large-scale structure and future sur- 
veys" in "Next Generation Wide-Field Multi-Object Spec- 
troscopy", ASP Conference Series vol. 280, ed. M.Brown & 
A.Dey (astro-ph/0301623) 
Eisenstein D.J., ct al., 2005, ApJ, 633, 560 
Eisenstein D..J., Sco H.-J., White M., 2007, ApJ, 664, 660 
Eisenstein D.J. et al., 2011, AJ submitted (arXiv: 1101.1529) 
Pukugita M., Ichikawa T., Gunn J.E., Doi M., Shimasaku K., 

Schneider D.P. 1996, AJ, 111, 1748 
Gilbank D.G., Gladders M.G., Yee H.K.C., Hsieh B.C., 2011, AJ, 
141, 94 

Glazebrook K., Blake C.A., 2005, ApJ, 631, 1 

Goliath M., Amanullah R., Astier P., Goobar A., Pain R., 2001, 

A&A, 380, 6 
Gunn J.E., et al., 1998, AJ, 116, 3040 
Gunn J.E., et al., 2006, AJ, 131, 2332 
Hamuy M., et al., 1996, AJ, 112, 2408 
Hamuy M., et al., 2006, PASP, 118, 2 
Hickcn M., et al., 2009, ApJ, 700, 331 

Hogg D.W., Finkbeiner D.P., Schlegel D.J., Gunn J.E., 2001, AJ, 

122, 2129 
Holtzman J.A., 1989, ApJS, 71, 1 
Holtzman J.A., et al. 2008, AJ, 136, 2306 
Hu W., Haiman Z., 2003, PRD, 68, 063004 
Hu W., Sugiyama N., 1996, ApJ, 471, 542 
Hutsi G., 2006, A&A, 449, 891 

Ivezic Z., ct al., 2004, Astronomische Nacliricliten, 325, 583 
Jarrctt T., Chester T., Cutri R., Schneider S., Strutskie M., 

Huchra J. P., 2000, AJ, 119, 2498 
Jha S., et al., 2006, AJ, 131, 527 
Jones D.H., et al., 2009, MNRAS, 399, 683 
Kaiser N., 1987, MNRAS, 227, 1 



20 Blake et al. 



Kazin E.A., ct al., 2010a, ApJ, 710, 1444 

Kazin E.A., Blanton M.R., Scoccimarro R., McBride C.K., 
Berlind A.A., 2010b, ApJ, 719, 1032 

Kazin E.A., Sanchez A.G., Blanton M.R., 2011, MNRAS submit- 
ted (arXiv: 1105.2037) 

Kessler R., et al., 2009, ApJS, 185, 32 

Komatsu E., ct al., 2009, ApJS, 180, 330 

Komatsu E., ct al., 2011, ApJS, 192, 18 

Kowalski M., et al., 2008, ApJ, 686, 749 

Krauss L.M., Turner M.S., 1996, GReGr, 27, 1137 

Landy S.D., Szalay A.S., 1993, ApJ, 412, 64 

Law N., et al., 2009, PASP, 121, 1395 

Lewis A., Challinor A., Lasenby A., 2000, ApJ, 538, 473 

Linder E.V., 2003, PRL, 90, 091301 

Lupton R., Gunn J.E., Ivezic Z., Knapp G.R., Kent S., 2001, in 
Astronomical Society of the Pacific Conference Series, Vol. 
238, Astronomical Data Analysis Software and Systems X, 
cd. F.R. Harnden, Jr., F.A. Primini, & H.E. Payne, 269-278 

Martin D., ct al., 2005, ApJ, 619, LI 
Matsubara T., 2008, PRD, 78, 083519 
McBridc C., ct al., in preparation 

Miller C.J., Nichol R.C., Batuski D.J., 2001, ApJ, 555, 68 
Ostriker J., Steinhardt P., 1995, Nature, 377, 600 
Padmanabhan N. et al., 2007, MNRAS, 378, 852 
Padmanabhan N., White M., 2008, PhRvD, 77, 3540 
Padmanabhan N., White M., 2009, PhRvD, 80, 63508 
Peebles P.J.E., Yu J.T., 1970, ApJ, 162, 815 
Percival W.J., et al., 2001, MNRAS, 327, 1297 
Percival W.J., ct al., 2007, ApJ, 657, 51 
Percival W.J., ct al., 2010, MNRAS, 401, 2148 
Pcrlmuttcr S., ct al., 1999, ApJ, 517, 565 

Pier J.R., Munn J. A., Hindsley R.B., Hennessy G.S., Kent S.M., 

Lupton R.H., Ivezic Z. 2003, AJ, 125, 1559 
Poole G., et al., in preparation 
Richards G., et al., 2002, AJ, 123, 2945 
Riess A.G., et al., 1998, AJ, 116, 1009 
Riess A.G., ct al., 1999, AJ, 117, 707 
Riess A.G., ct al., 2004, ApJ, 607, 665 
Riess A.G., ct al., 2007, ApJ, 659, 98 
Rimes C.D., Hamilton A.J.S., 2005, MNRAS, 360, 82 
Ross A., et al., 2011, MNRAS submitted (arXiv: 1105.2320) 
Sanchez A.G., Baugh C.M., Angulo R., 2008, MNRAS, 390, 1470 
Sanchez A.G., Crocce M., Cabre A., Baugh G.M., Gaatanaga E., 

2009, MNRAS, 400, 1643 
Saunders ct al., 2004, SPIE, 5492, 389 
Sec H.-J., Eisenstein D.J., 2003, ApJ, 598, 720 
Seo H.-J., Siegel E.R., Eisenstein D.J., White M., 2008, ApJ, 686, 

13 

Sharp et al., 2006, SPIE, 6269, 14 

Shoji M., Jeong D., Komatsu E., 2009, ApJ, 693, 1404 

Smith J. A., et al., 2002, AJ, 123, 2121 

Smith R.E., et al., 2003, MNRAS, 341, 1311 

Smith R.E., Scoccimarro R., Sheth R., 2008, PhRvD, 77, 43525 

Stoughton G., et al., 2002, AJ, 123, 485 

Sunyaev R.A., Zeldovitch, Y.B., 1970, Ap&SS, 7, 3 

Takahashi R., ct al., 2011, ApJ, 726, 7 

Taruya A., Nishimichi T., Saito S., 2010, PRD, 82, 3522 

Tegmark M., 1997, PRL, 79, 20 

Thomas S., Abdalla F., Lahav O., 2011, PhRvL, 106, 1301 
Tucker D.L., et al., 2006, Astronomische Nachrichten, 327, 821 
White M., et al., 2011, ApJ, 728, 126 
Wood-Vasey W.M., et al., 2007, ApJ, 666, 694 
York D.G. et al., 2000, AJ, 120, 1579 



